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ABSTRACT 

Context. We present a detailed NLTE analysis of 39 Mn i lines in the solar spectrum. The influence of NLTE effects on the line 
formation and element abundance is investigated. 

Aims. Our goal is the derivation of solar log gfe values for manganese lines, which will later be used in differential abundance analy- 
ses of metal-poor stars. 

Methods. The method of spectrum synthesis is employed, which is based on a solar model atmosphere with initially specified element 
abundances. A manganese abundance of log fijvm.o = 5.47 dex is used with the theoretical line-blanketed model atmosphere. Statistical 
equilibrium calculations are carried out for the model atom, which comprises 245 and 213 levels for Mm and Mnii, respectively. 
Photoionization cross-sections are assumed hydrogenic. 

Results. For line synthesis van der Waals broadening is calculated according to Anstee & O'Mara's formalism. It is shown that hyper- 
fine structure of the Mn lines also has strong broadening effects, and that manganese is prone to NLTE effects in the solar atmosphere. 
The nature of the NLTE effects and the validity of the LTE approach are discussed in detail. The role of photoionization and collisional 
interaction is investigated. 

Conclusions. Maximum NLTE corrections of +0.1 dex with respect to LTE profiles are found. We propose an absolute solar abun- 
dance of 5.36 + 0.1 dex. The main source of errors in the abundance calculations is uncertain oscillator strengths. 
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1. Introduction 

Manganese belongs to the iron-group elements, which have very 
complex but interesting atomic structure. Since lines in the syn- 
thetic spectrum of a star are based on the adopted model atom 
and the transitions included, the analysis becomes more and 
more complicated for elements with increasing number of lev- 
els and electrons. Construction of the model atom, usually done 
within a few days, becomes a matter of several months. It is nec- 
essary not only to collect all existing data and compare them, 
but also to assess their plausibility and necessity for the aims of 
current research. 

These difficulties are, however, compensated by a beautiful 
pattern of lines in the spectrum, which results from fine and hy- 
perfine transitions in ions of such elements. The Mm lines are 
split into several components; this is a result of interaction of the 
electromagnetic field produced by the electrons with the mag- 
netic momentum of the nucleus. The large hyperfine structure 
(HFS) of the lines overwhelms all other sources of line broaden- 
ing. The width of a single Doppler line at 8400 A in the solar at- 
mosphere is ~ 28 mA versus ~ 40 mA of a HFS broadened line. 
Without question, it has to be accurately taken into account, oth- 
erwise serious errors appear in the calculation of element abun- 
da nces (A bt [T952l Proc haska & McWiniam |2000l del Peloso et 
al. 12005'; Vitas & Vince!2003]l. 

It is interesting to study the manganese abundance in stars of 
different populations. In metal-poor stars the manganese abun- 
dance was found to be correlated with metallicity. A progres- 
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sive deficiency of Mn relative to iron was first reported by 
Heifer et al. ( 1 1 9591 ) and later confirmed in analyses of other au- 
thors (Wheeler et al. T9 891 Gr atton ,1989; Nissen et al. l2000l 
Prochaska & Mc William 120001) . right up to the most recent in- 
vestigation of Sobeck et al. (2005). This element is not produced 
in quasi-static nuclear burning processes; instead, it is accumu- 
lated from supernova explosions. Therefore, the variation of the 
abundance pattern of manganese with metallicity may indicate 
when and which stars (SN la, SN II), and which production pro- 
cesses of heavy elements started to contribute to the chemical 
enrichment of the Galaxy. 

Aside from global implications for studies of stellar nucle- 
osynthesis, manganese is an indicator of NLTE conditions in a 
stellar atmosphere. Mn does not contribute to stellar opacity, as 
compared to iron, and it does not significantly add to the pool of 
free electrons. More than 95% of manganese is ionized. Thus, 
Mn I can be regarded as a trace ion in a stellar atmosphere and 
studied under the NLTE assumption, which is necessary to de- 
scribe lines formed high in the photosphere (Mihalas & Athay 
II973I I. Departures from LTE for the resonance Une at 5394 A 
in the Sun have been already found by Vitas & Vince ( 120051 1. 
although the nature of the NLTE effects was not investigated. 
Nissen et al. (2000) proposed that Mn abundances might be sub- 
ject to NLTE effects due to overionization of Mm caused by a 
strong radiation field. 

The goal of the cuiTent analysis is to investigate in detail de- 
partures from LTE in Mn i and assess the plausibility of previous 
element abundance estimates for the Sun. There is an unresolved 
inconsistency between the abundance of Mn in the solar atmo- 
sphere (5.39 + 0.03) (Booth et al. |1984bl l and in CI Chondrites 
(5.50 + 0.03) (Lodders 2003j. It is not yet clear whether to at- 
tribute this to an incorrect allowance for HFS and NLTE effects 
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in the calculation of the solar Mn abundance. The fit of photo- 
spheric and meteoritic abundances of the reference element Si is 
also very important in this context. 

In Sect. 2 we give a short account of the method used and 
the full atomic model of Mm including the choice of interac- 
tion processes. For the purpose of identifying the important lev- 
els and transitions, this section also presents two reduced mod- 
els. Section 3 explores the resulting departure coefficients, line 
source functions, and the nature of the deviations from LTE. It 
also introduces their variation with photoionization and collision 
cross-sections. In Sect. 4 we dicuss the synthesis of the Mn i lines 
with due reference to HFS and NLTE results. The resulting solar 
manganese abundance is discussed. 

2. NLTE calculations 

2.1. Method 

In our current atomic model the following processes are taken 
into account: photon absorption in line transitions, photoioniza- 
tion, excitation and ionization by collisions with free electrons, 
and neutral hydrogen atoms. All processes include their reverse 
reactions; in particular, bound-bound radiative transitions are as- 
sumed to follow complete frequency redistribution. It is known 
that under certain circumstances other interaction processes can 
also play a role, e.g. autoionization transitions, charge-transfer 
reactions or di-electronic recombination. No information about 
these processes is available for manganese, and in particular, au- 
toionization resonances could well contribute to the depopula- 
tion of some Mm levels. 

In this work NLTE atomic populations were computed with 
a revised version of the program DETAIL (Butler & Giddings 
11985 ), which solves the radiative transfer and statistical equi- 
librium equations by the method of accelerated lambda itera- 
tion. We use a theoretical line-blanketed model of the solar pho- 
tosphere MAFA GS-OD F, calculated with the MAFAGS code 
(Fuhrmann et al. II997I I. This model atmosphere uses opacity 
distribution functions (ODF) taken from Kurucz (1992). The re- 
sulting atmospheric stratifications of temperature and pressure 
are similar to those given by other comparable models (see com- 
parison in Grupp 2004, Fig. 15). As is the case with all other 
line-blanketed atmospheric models of this type so far, we have 
not attempted to model the solar chromosphere. 

2.2. Models of the Mn i atom 
2.2.1. The reference model 

In our reference model the manganese atom is constructed with 
245 levels for Mm and 213 levels for Mnii, respectively. The 
system is closed with the ground state of Mn m. The energies for 
these levels are taken from Sugar & Corliss ( 119851 ). We include 
all levels with n < 15 for Mm, n < 8 for Mmi, and with ener- 
gies of 0.03 and 1.11 eV below the respective ionization limits. 
Such a relatively complete model should provide a very close 
coupling of the upper atomic levels to the next ion ground state. 
The number of lines treated in NLTE is 1261 for Mm and 1548 
for Mn 11. 

Wavelengths and oscillator strengths are all taken from 
Kurucz's database (Kurucz & Bell I19951 I. We note here that 
most of the lines do not require data that are more precise (in 
fact it will turn out that the available experimental f-values are 
not better than those provided by calculations). In a few cases 
of strong line blending in the terrestrial UV, we have tried to 



make sure that the calculated wavelengths are at the observed 
positions. Hyperfine splitting of the lines is not included in the 
statistical equilibrium calculations, since there is no reason to 
believe that the relative populations of these levels deviate from 
thermal. Moreover, for the uppermost terms above 7 eV the fine 
structure is not maintained. These terms are represented by a 
single level with a weighted mean of statistical weights and ion- 
ization frequencies of their fine structure levels. A complete gro- 
trian diagram for Mm is available online. The majority of radia- 
tive transitions in Mn i occurs between low excitation metastable 
and high excitation levels. Also, transitions starting from the 
metastable levels produce most of the lines, which we use in the 
spectrum analysis. Quite different from Fe i, the intercombina- 
tion lines do not seem to form a tight coupling of the multiplets. 
Except for the calculation of the ionization equilibrium, the mul- 
tiplets could therefore be calculated for themselves ignoring any 
multiplet interaction. 

For bound-free radiative transitions we have to use 
hydrogen-like photoionization cross-sections (Mihalas 119781 ) 
because calculated data from the Opacity Project are not yet 
available. In our current analysis this may be the most uncer- 
tain representation. We conclude this in analogy to the Fe i atom, 
for which we found that the calculated photoionization cross- 
sections of Bautista ( 1 19971 1 are orders of magnitude larger than 
hydrogenic approximations (Gehren et al. 12001 at . 

Rates of transitions due to inelastic collisions with H i atoms 
are calculated according to Drawin's formula ( 119681 1 19691 ) in 
the version of Steenbock & Holweger ( II984I ). Drawin's cross- 
sections are usually multiplied by a scaling factor 5h, which 
often takes values much smaller than unity (see the review by 
Asplund 2005 ). Thus, we decided to scale the rates of bound- 
bound and bound-free transitions by a factor of 0.05 in our refer- 
ence model. For allowed transitions, the cross-sections for colli- 
sional excitation with electrons are calculated from the formula 
of van Regemorter (119621 ). Forbidden transitions due to colli- 
sions with electrons are computed from the formula of Allen 
(fT973]) . 

2.2.2. Reduced models 

Our procedure of reducing the complexity of the model atom 
lies in a consecutive exclusion of those atomic levels that are 
suspected to cause the secondary NLTE effects, such as photon 
pumping or radiative recombination. These levels are identified 
in Sect. 3.1. The purpose is to find a model that represents the 
idealized case of the manganese atom, where deviations from 
LTE in level populations are exclusively due to one primary pro- 
cess. We will also be able to establish a minimum model that 
gives similar results as the reference one. As a by-product the 
time of computation could be significantly decreased. 

As high-excitation levels provide an effective coupling to the 
continuum by means of collisions, it is tempting to exclude them 
in order to minimize the collisional interaction of Mn i and Mn ii. 
A first reduced model of Mm thus includes all 145 levels with 
an excitation potential of less than 6.42 eV {Eion - 7.43 eV). 
The Mm levels with n - 6 ... 15 and all levels of Mnii, except 
the ground state, were excluded from the calculations. 

Since the increase to 1 eV of the energy gap between both 
ions is not enough to ensure the dominance of radiative pro- 
cesses, and it does not lead to perceptible changes in the distri- 
bution of level populations (see next section), a second reduced 
model was calculated with only 65 levels (n < 5) for Mm, and 
an excitation potential of the uppermost level of 5.23 eV. Both 
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types of reduced models in practice eliminate the complete dou- 
blet system except for the metastable states. 

Whereas the above reduced model atoms are defined with re- 
spect to the excitation energies of the contributing levels, other 
models may exclude single multiplet systems, such as the dou- 
blet or the octet system. Note that statistical equilibrium requires 
a valid representation of such levels in order to calculate the ion- 
ization equilibrium. In the case of the doublet system it seems 
that it can be completely ignored because it contains only a few 
levels and all of them are highly excited. Also, there are no use- 
ful doublet lines in the solar spectrum. This is different with 
the octet system, which is also only loosely coupled to quartets 
and sextets. However, a significant fraction of the Mm atoms is 
found here, in addition to some interesting solar lines. 

Although we find that the number of levels (and lines) in 
both Mn I and Mn ii can be significantly reduced without chang- 
ing the formation of the most important solar manganese lines, 
we have calculated all final results with the complete reference 
model from Sect. 12.2.11 The reduced models are therefore only 
used for direct investigation of critical processes. 



3. Results 

3.1. NLTE level populations 

Departure coefficients bi = nj^^™/n['™ for Mm levels, calcu- 
lated under different model assumptions, are presented in Figs.[T] 
and|2l where Fig.[T^ shows the deviations from LTE for the ref- 
erence atom. Each term is represented by one level because of 
the close coupling among the fine-structure levels. In order not 
to overload the figure, we have only entered a few selected lev- 
els which are typical for their depth dependence. As an example, 
the relative population of the ground state a ^S^^j is very simi- 
lar to that of the other metastable terms, a^D, a'^D, and a'*G. 
The thick curves refer to levels of relatively low excitation; they 
are confronted with a number of highly-excited levels that show 
a very different behaviour These levels deviate from the low- 
excitation levels following a trend that predicts more thermal or 
even superthermal populations the higher the excitation energies 
are. 

In the deepest layers, b, - 1 for all levels (see Fig. [T^). 
Deviations from LTE develop between logTsooo ~ 0.1 ... 1. At 
these depths the atmosphere becomes optically thin below the 
ionization thresholds of some low-lying levels, such as a'^D 
(A < 2729 A), z«P° (A < 2844 A), y^P° (A < 2884 A), and 
a'^H (A < 3876 A). The detailed balance is broken because pro- 
cesses of radiative ionization start dominating over recombina- 
tions, which depend on the local temperature. Depopulation of 
these levels is followed by a depopulation of the Mm ground 
state. This is powered by strong collisional processes, which 
couple fl*Sg^2 '^o the excited levels (essentially the metastable 
ones) and maintain the relative balance among them out to 
logrgooo ~ -2.8. 

Radiation in the stronger lines also contributes to the re- 
distribution in the atomic level populations. Line pumping 
takes action when the optical depth drops below unity in the 
wings of lines that belong to multiplets 21, 22 and 23 ( 
A 4257,4761,4739,4451 A). A strong non-local UV radiation 
field, Jy > By{Ts), increases the photoexcitation rates, which are 
not compensated by the inverse rates of de-excitation. Hence, at 
the depths of line formation, below log T5000 ~ -1.5, populations 
of the upper levels for these transitions (y^P", z'^D", z'^F") are 



amplified relative to the ground state and adjacent low excitation 
levels. The line source functions S ij are larger than By{T^). 

The photon suction (Bruls et al. 119921 effectively operates 
in the strong lines of multiplets 27 (transitions from e^S to 
z^P°) and 16 (e^^S to z^P°). This process can be understood 
as a successive radiative de-excitation of an atom through tran- 
sitions with high probability in atmospheric layers where the 
mean intensity Jy falls below By(T^), resulting in a net down- 
ward rate. This leads to a perceptible depopulation of the upper 
levels e and e ^S, although the corresponding overpopulation 
of the lower levels is rather small. 

Levels with an excitation energy above 6.5 eV are more 
closely linked to the continuum by means of collisions. For the 
majority of these high -excitation Mn i levels the processes of re- 
combination will dominate instead, which is evident from their 
negative net ionization rates. The reason is that at large depths 
Jy in the IR spectrum is smaller than By(Te). Lack of ionizing 
photons will result in an increased net recombination. Hence the 
levels will be slightly overpopulated, most of them remaining in 
detailed balance with a ^83. 



3.2. Nature of the NLTE effects 

It is clear from the analysis given in the previous section that 
there are several processes that take part in establishing the oc- 
cupation numbers of Mm and Mnn levels. To identify the ma- 
jor cause of the non-equilibrium populations, we perform NLTE 
calculations applying various scaling factors to photoionization 
and collision cross-sections. In addition, we investigate how the 
completeness of the model atom affects the results. 



3.2.1. Standard models 

The reference atomic model of Mn i, as described in Sect. 12.2.11 
was constructed with hydrogenic photoionization cross-sections 
and cross-sections for inelastic collisions with H i (both for b- 
b and b-f transitions) based on Drawin's formula, but scaled by 
a factor Sh = 0.05. Comparison of the collision rates C,y(H) 
and C,y(e) shows that in the solar atmosphere for Mm rates of 
collisions with electrons are on average greater by one order of 
magnitude. In order to test the influence of inelastic collisions 
with hydrogen, we carried out additional calculations for scaling 
factors from 5h = to 5 (Fig. [T^-c). Naturally, with increas- 
ing rates of collisional excitations and ionizations by H i atoms 
departures get closer to unity, although perfect thermalization is 
not achieved with values for as low as 5. 

The influence of photoionization was studied likewise, i.e., 
by means of performing NLTE calculations for various scaling 
factors to the cross-sections. For test purposes we have chosen to 
scale the Mn i cross-sections cr,, by factors Sp = 0, 300 and 5000. 
A plot of the resulting departure coefficients is presented in Fig. 
[TJl-f. Absence of photoionization (Fig.[TJl) does not restore LTE, 
as we could expect if the major and only NLTE mechanism were 
the overionization. Moreover, the results for zero photoioniza- 
tion cross-sections are virtually indistinguishable from those of 
the reference model (Sp - I), while Sp - 300 leads to slightly 
stronger deviations of bj from unity (Fig. [1^). When cr^ are in- 
creased by a factor of 5000 (Fig. [if) the overionization is ampli- 
fied and all Mn i levels are noticeably affected. 
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Fig. 1. Departure coefficients b, of selected Mm levels for our standard model atom as a function of optical depth. The curves are 
labeled on the right, (a): The reference model atom with a total of 459 levels. Hydrogen collision rates are scaled by Sh - 0.05. 
The ground state a ^S^^j representative of most of the metastable levels such as a ^D, a '^D, and a '^G. The odd levels z **P° and 
z^P° are also very similar. Note a general trend of the departure coefficients increasing with excitation energy, (b): Same as (a), but 
hydrogen collisions set equal to zero. The departure coefficients are very much the same as in the reference model, (c): same as 
(a), but Sh - 5.0. (d): same as (a), but photoionization rates set equal to zero.(e); same as (a), but hydrogenic photoionization rates 
scaled by a factor 300. (f): same as (a), but photoionization rates scaled by a factor 5000. 



3.2.2. Reduced models 

We consider reduced model atoms that are discussed in detail 
in this section. In Fig.|2]we show departure coefficients for two 
reduced models, one with 145 levels of Mm and an excitation 
potential of the uppermost level of 6.42 eV (/sion - 7.43 eV), 



and one with 65 Mn i levels limited to 5.23 eV. For both models 
only the Mn ii ground state was considered. Note that due to the 
model reduction some of the high-excitation levels in Fig. [T] are 
no longer present. 
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Fig. 2. Departure coefficients b, of selected Mm levels for reduced model atoms as a function of optical depth. The curves are 
labeled on the right, (a): The reduced model atom with a total of 146 levels including the Mnii ground state. Hydrogen collision 
rates are scaled by Sh - 0.05 as in the reference atom. Note that the uppermost levels of Fig. [1] are no longer represented by the 
reduced model, (b): Same as (a), but photoionization rates set equal to zero, (c): The reduced model atom with a total of 66 levels 
including the Mnii ground state. Hydrogen collision rates are scaled by Sh = 0.05 as in the reference atom, (d): same as (c), but 
photoionization rates set equal to zero, (e); same as (c), but electron collision rates set equal to zero, (f): same as (d), but electron 
colUsion rates set equal to zero. 



Fig. and b display the departure coefficients for the 146 
level atom. As a result of a reduced recombination efficiency the 
remaining high-excitation levels (not shown here) are slightly 
less overpopulated, when compared to the NLTE populations 
calculated with the full atomic model. Fig. |2j5 shows calcu- 



lated with zero photoionization. Most of the low-excitation lev- 
els are still underpopulated due to pumping processes combined 
with an efficient collisional coupling to the Mn ii ground state. It 
is only for the 66 level model atom that the domination of pho- 
toionization can be seen. Fig. |2j; and d demonstrate the loss of 
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the hydrogen collision efficiency in coupling to the Mn n ground 
state. Here, removal of photoionization with Sp = drives the 
lower level populations towards LTE, at least in the line forma- 
tion region. In the absence of photoionization the photon suc- 
tion along the lines now results in an overpopulation of low- 
excitation Mn I levels that is balanced by net collisional ioniza- 
tion. Neglecting electron collisions (Fig.|2^ and f) results in an 
even weaker coupling of the levels. 

3.2.3. Discussion 

The behaviour of Mn level populations, strongly deviating from 
LTE, can not be explained simply in terms of photoionization or 
collisions, as was succesfully done for other elements submit- 
ted to an NLTE analysis. Our results show that hydrogenic b-f 
radiative cross-sections for the Mn i levels of 3 to 4 eV excita- 
tion are too small to ensure the dominance of photoionization, 
which is the main NLTE effect for many metal-group elements. 
These levels usually dominate the ionization balance in such ions 
because in the atmospheres of solar-type stars their photoion- 
ization edges lie in regions where the radiation field is intense 
enough and the radiative rates prevail over collisional rates. For 
instance, the Opacity Project calculations (Bautista '1997') give 
evidence for neutral iron cross-sections, which are 10 to 1000 
times larger than the hydrogenic ones. This could explain the 
nearly identical behaviour of Mn departure coefficents for 5p = 
and 5p = 1 . Kramers' approximation for Mn i gives a lower limit 
for the cross-sections, which is comparable in its influence on 
populations to the absence of photoionization. Thus, a very large 
enhancement to o-y is required in order to force the dominance 
of this process in establishing the occupation numbers of Mn i 
and Mnii levels. 

The complexity of the Mm atomic system with its large 
number of levels produces some additional complications for 
NLTE calculations. There are many doubly excited (DE) config- 
urations, which, in contrast to those of single excitation (SE), can 
not be treated in hydrogenic approximation. The use of Kramers' 
formula with principal quantum number n for all atomic levels 
is an obvious simplification: DE and SE configurations of neigh- 
bouring excitation energies but different n will have strongly 
differing photoionization cross-sections, which scale with in- 
verse fifth power of the principal quantum number Furthermore, 
the b-f collision cross-sections, calculated using the radiative 
b-f cross-sections, will differ for these neighbouring levels by 
some orders of magnitude. The latter is important for the high- 
excitation levels whose b-f collision rates compensate radiative 
b-f rates from the low-excitation levels and in this way establish 
the statistical equilibrium of the model atom. As seen on Fig.|2}; 
and d exclusion of the upper levels with ~ 1 - 2 eV energy sep- 
aration from the continuum, leads to a reduced recombination 
efliciency, and the Mn i atom could be identified as a member of 
the photoionization-dominated category. 

One of the possible approaches is to use the effective princi- 
pal quantum numbers of the levels. We have performed NLTE 
calculations for the photoionization cross-sections computed 
from Kramers' formula with Wetf ■ The switch from principal to 
effective quantum numbers did not change the results signifi- 
cantly except that the uppermost levels are now coupled more 
tightly to the Mnii ground state at the line formation depths. 
Also, the whole atomic system became more sensitive to the 
enhancement of photoionization cross-sections, but levels that 
produce the lines of interest in our research were not affected. 
This solution may be more appropriate for very complex atomic 
systems, such as Mn, but at the same time it emphasizes the inad- 



equacy of the hydrogenic approximation. This is, however, noth- 
ing new; it was also shown by Nahar & Pradhan (2004) that pho- 
toionization of excited levels of a non-hydrogenic ion may not 
be treated in hydrogenic approximation. As there is currently no 
alternative to the latter, we decided to preserve the hydrogenic 
photoionization cross-sections. 

The incompleteness of the available data for the transitions 
in the Mn model atom presents another problem. There are 
more than 500 experimental levels reported for Mm, 180 being 
above the lowest ionization limit (Johansson & Cowley 1988l l. 
Although the system of atomic levels seems to be complete, the 
case is somewhat different for the allowed transitions between 
them. In particular, there is no data in Kurucz's database (Kurucz 
& Bell |1995l l, which we use in our analysis, for transitions that 
link the uppermost levels and levels of intermediate energies. 
Unfortunately, it is exactly these transitions that should couple 
the bulk of the overpopulated levels to the remaining underpop- 
ulated levels in our model. 

In order to check the results obtained with the DETAIL code, 
we performed statistical equilibrium calculations for Mn using 
the independent N0NLTE3 code, based on the method of com- 
plete linearization (Sakhibuflin 1983). The NONLTE3 atomic 
level populations are in close agreement with those from the 
DETAIL code. Hence, we do not associate the behaviour of the 
Mn atomic system with program errors. We suggest that the Mn i 
ion represents a complex mixed-domination case, where inelas- 
tic collisions with electrons and hydrogen atoms counteract pho- 
toionization and decouple the atomic system. 

4. Synthesis of Mn i lines 

Spectrum synthesis is employed to determine the abundance of 
Mn in the solar atmosphere. We use our line analysis code SIU, 
which generates synthetic NLTE line profiles for computed level 
departure coefficients ^/(t). For all elements, except Mn, we as- 
sume LTE. 

The solar spectrum is calculated using the MAFAGS-ODF 
solar model atmosphere with Teff = 5780 K, an initial man- 
ganese abundancqjof log fiMn.o - 5.47 dex, and a constant mi- 
croturbulence velocity ft = 0.9 km s"'. The profiles are broad- 
ened by a rotation velocity Vrot = 1 .8 km s"' , and by a macrotur- 
bulence velocity Vmac = 2.5 ... 4 km s ' . The observed spectrum 
is taken from the Kitt Peak Solar Flux Atlas (Kurucz et al. |1984| l. 
For the abundance analysis we originally selected 39 Mn i lines 
that satisfy the following conditions: they are relatively free from 
blends, and laboratory data for hyperfine splitting and oscillator 
strengths are available. Relative intensities of the HFS compo- 
nents are calculated according to the tables of White & Eliason 
C1933'). Magnetic dipole splitting constants, A{J), and electric 
quadrupole splitting constants, B(J), of the corresponding levels 
are given in the online table. Van der Waals' damping constants 
logCe are computed according to Anstee & O'Mara's ( 1995l l 
formalism; a correction of A log Ce = -0. 1 is applied in order 
to fit the wings of strong lines. We use oscillator strengths from 
various sources giving preference to measured values. All in- 
put parameters necessary to perform spectrum synthesis for the 
lines of interest are given in Table[T] Line profiles are computed 
under both LTE and NLTE assumptions; they are fitted to the 
observed profiles by means of manganese abundance variations. 
The logarithmic abundance differences of the fits with respect 
to the initial manganese abundance, A log e, are reproduced in 

' Note that the atmospheric model of the Sun does not depend on the 
manganese abundance. 
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Table 1. Parameters of lines used for abundance calculation. An asterisk (*) in the wavelength entry refers to the lines with consistent 
ocillator strengths from different sources (see text). No error estimate is given for Kurucz \og gf values. 



No. 


A 


Mult. NhfS ^low 


Lower 


Upper 


Wa log gf Error 


A lege 


AX 


log g/e Source^ 




[A] 


[eV] 


level 


level 


[mA] 









136. 


-0.070 


0.1 


-0.15 


0.0 


5.250 


101. 


-0.446 


0.1 


-0.05 


0.0 


4.974 


70. 


-0.950 


0.1 


-0.06 


0.06 


4.460 


62. 


-0.700 


0.18 


0.11 


0.095 


4.880 


71.3 


-0.288 


0.04 


-0.21 


0.07 


4.972 


93. 


0.278 


0.04 


-0.3 


0.0 


5.448 


53.5 


-0.490 


0.04 


-0.165 


0.08 


4.815 


73. 


-0.246 




0.00 


0.05 


5.224 


50. 


-0.555 


0.1 


-0.185 


0.06 


4.730 


53.5 


-0.444 


0.04 


-0.23 


0.07 


4.796 


56. 


-0.522 


0.04 


-0.16 


0.065 


4.788 


57. 


-0.343 


0.04 


-0.37 


0.06 


4.757 


59. 


-0.345 


0.04 


-0.37 


0.05 


4.755 


12.8 


-1.670 


0.1 


-0.13 


0.08 


3.670 


72. 


-0.340 


0.04 


-0.35 


0.06 


4.780 


62. 


-0.490 


0.04 


-0.24 


0.07 


4.740 


146. 


-0.086 


0.04 


-0.16 


-0.04 


5.224 


73. 


-0.138 


0.04 


-0.18 


0.07 


5.152 


108. 


0.425 


0.04 


-0.36 


0.0 


5.535 


81. 


-0.080 


0.1 


-0.17 


0.05 


5.220 


98.5 


0.100 


0.1 


-0.23 


0.0 


5.340 


148. 


0.042 


0.04 


-0.32 


0.0 


5.192 


149. 


0.144 


0.04 


-0.32 


-0.03 


5.294 


13.7 


-1.630 


0.1 


-0.12 


0.08 


3.720 


24.2 


-1.140 


0.1 


-0.05 


0.08 


4.280 


41.5 


-0.763 


0.04 


-0.2 


0.08 


4.507 


79.5 


-3.503 


0.1 


-0.115 


0.105 


1.852 


53. 


-1.743 


0.1 


-0.125 


0.085 


3.602 


85. 


-1.462 


0.1 


-0.06 


0.09 


3.948 


48. 


-3.795 


0.1 


-0.12 


0.11 


1.555 


57.5 


-1.702 




-0.04 


0.085 


3.728 


44. 


-1.847 


0.1 


-0.02 


0.085 


3.603 


36. 


-2.017 




-0.02 


0.085 


3.433 


87. 


-0.251 


0.1 


-0.285 


-0.015 


4.934 


97.8 


-0.216 




-0.145 


-0.05 


5.109 


96.8 


0.034 


0.1 


-0.27 


-0.07 


5.234 


11. 


-0.477 




-0.058 


0.1 


4.935 


15.5 


-0.328 




-0.057 


0.1 


5.085 


26. 


-0.055 




-0.09 


0.09 


5.325 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 



4055.513* 

4058.911* 

4070.264 

4257.641 

4436.342 

4451.581 

4453.001 

4455.288 

4457.010* 

4470.143 

4490.067 

4498.901 

4502.220 

4671.667* 

4709.705 

4739.088 

4754.021* 

4761.508 

4762.358 

4765.851 

4766.413 

4783.389* 

4823.460* 

5004.894 

5117.913* 

5255.287* 

5394.619* 

5407.331* 

5420.265* 

5432.512* 

5470.560 

5516.697* 

5537.692 

6013.465 

6016.586 

6021.727 

8700.877 

8703.559 

8740.705 



5 
5 
5 
23 
22 
22 
22 
28 
28 
22 
22 
22 
22 
21 
21 
21 
16 
21 
21 
21 
21 
16 
16 
20 
32 
32 
1 
4 
4 
1 
4 
4 
4 
27 
27 
27 
49 
49 
49 



4 
3 
3 
3 
3 
3 
2 
3 
3 
2 
2 
2 
2 
5 
4 
4 
5 
4 
5 
3 
4 
5 
6 
4 
3 
6 
6 
10 
9 
5 



5 
6 
6 
6 
10 
11 
13 



2.13 
2.17 
2.19 
2.94 
2.91 
2.88 
2.93 
3.07 
3.06 
2.93 
2.94 
2.93 
2.91 
2.88 
2.88 
2.93 
2.27 
2.94 
2.88 
2.93 
2.91 
2.29 
2.31 
2.91 
3.12 
3.12 


2.13 
2.13 


2.15 
2.17 
2.18 
3.06 
3.06 
3.06 
4.41 
4.41 
4.42 
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,4 

-"3/2 
^%2 

7/2 
3/2 
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am 
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fl'*G,,/2 
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3/2 
1/2 



^^3/2 
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y '^5/2 
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10 
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1/2 
. "2 

^ "^3/2 

^ '-'7/2 

^ '-'1/2 

" '-'312 

e U3/2 

^ ^-^3/2 

^ '-'3/2 

^ ^^5/2 

^ '-'7/2 

^ ^5/2 

'o 

7/2 
3/2 

7/2 

'o 

3/2 
7 4po 
i. 1 g/2 

^ ^^^5/2 
i. 1 7/2 

e*S 

e O7/2 
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z^F° 

i. 1 7/2 
,4po 
^ ^3/2 
,4po 
1 9/2 
8po 
<. 1:7/2 
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References: (1) Booth et al. |1984al (2) Greenlee & Whaling [T979l (3) Blackwell & Collins [T9721 (4) Woodgate [T966l (5) Kurucz [T988] 



Table [U (Col. 11). In this paper we refer to this parameter as 
abundance correction. The difference in abundances required to 
fit LTE and NLTE profiles is referred to as the NLTE correction 
(AX = log e^^LTE _ ^TE). gjygjj (-qJ ^2 of Table[I]for 

each line. 



4.1. Line profiles 

The computed profiles for selected lines are given in Fig. [3] to- 
gether with the solar spectrum. For comparison we show profiles 
under LTE and NLTE conditions, with and without HFS. The ex- 
amples on the left and on the right panel are representatives of 



two groups of lines that we were able to distinguish in the anal- 
ysis. 

The first group contains those lines, mainly weak with W,\ < 
80 mA that are formed deep in the photosphere, in the layers 
with a strong gradient of the local temperature Tg. The forma- 
tion of lines under NLTE takes place at higher when com- 
pared to the case of LTE. As bi < bj over the entire line forma- 
tion depth, the line source function is larger than its LTE value, 
S ij > By{Tfi). Consequently, NLTE corrections AX will be pos- 
itive due to decreased absorption over the entire line profile. In 
this case, both LTE and NLTE profiles will fit the observed lines, 
provided a certain abundance correction is performed. Of course. 
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Fig. 3. Selected observed solar Mm profiles (filled circles). Synthetic NLTE and LTE profiles with and without HFS are labeled 
correspondingly. NLTE results are based on the reference atomic model. Wavelength positions and relative line strengths of the 
HFS components are indicated. The line abundances and van der Waals' damping constants logCe are -0.02,-31.41 (5516 A); 
-0.32,-30.78 (4783 A); -0.13,-30.73 (4671 A); -0.27,-30.64 (6021 A); -0.11,-31.40 (5394 A); and -0.16,-30.80 (4754 A), 
where the line abundance parameter refers to the abundance difference with respect to logSMn,© = 5.47 dex. 



abundances obtained from the LTE profile fitting will be erro- 
neous. 

Stronger lines are formed within very large depth ranges. At 
the depths of line core formation the upper levels are depopu- 
lated more efficiently than the lower ones, h\ > bj. The source 
function drops below its LTE value, S ij < ByiJe), which is 
driven by photon escape in the line wings. In the lower atmo- 
sphere, behaviour of the source function is inverted, i.e., bi < bj. 
Accordingly, for the strongest lines (e.g., AA 4055, 4783, 6021 A) 
we observe an amplified absorption in the core and a decreased 
absorption in the wings relative to LTE. It is important that pro- 
files computed under the LTE approach can not fit the observed 
lines due to their different profile shapes. 

It is noteworthy that we succeeded to fit all lines in the list 
with a single value of microturbulence velocity Although it is 
well-known that this parameter approximates a depth-dependent 
velocity field in the solar atmosphere, we decided to preserve a 
constant = 0.9 km s"' . Certainly, this assumption does not al- 
ways give the best result: lines of intermediate strength require 



different values of the microturbulent velocity to fit the core and 
the wings simultaneously. We have not tried to adjust this value. 
Our analysis of the solar lines is mainly aimed at derivation 
of atomic and atmospheric data to be used in metal-poor stars, 
where only a single depth-independent value of can be recov- 
ered. Other asymmetries in the observed lines, such as red wing 
deficit or core skewness (as seen on all lines in the right panel), 
are a reflection of atmospheric inhomogeneities, which we can 
not take into account with our static plane-parallel models. 

4.2. Solar abundance of Mn and oscillator strengths 

The method of spectrum synthesis yields the product of the os- 
cillator strength for a transition and the abundance of the ele- 
ment, log gfs. This parameter, as well as NLTE correction AX, 
and abundance correction A logs, are given in Cols. 13, 12, 
and 1 1 of Table [1] respectively. Figure |4^ shows NLTE and 
LTE abundances for all lines as a function of their equivalent 
widths. It is seen that the absolute abundances determined as- 
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Table 2. Weighted NLTE manganese abundances based on os- 
cillator strengths from different sources. 



^lines 


Source of 


Method 


accuracy 


logs 


cr 




ioggf" 




% 


[dex] 




17 


1 


beam foil 


25 


5.35 


0.07 


20 


2 


beam foil 


50 


5.46 


0.2 


19 


3 


laser excitation 




5.30 


0.08 


25 


4 


total absorption 


10 


5.23 


0.10 


39 


5 


calculations 




5.37 


0.08 


39 


all 






5.28 


0.11 



References: (1) Green lee & Whaling [T9791 (2 ) Woodgate 1966.; 
(3) Becker et al. [T980l (4) Booth et al. |1984al (5) Kurucz [19881 



suming LTE are lower than those derived with NLTE level popu- 
lations. We find the weig htecQ mean NLTE and LTE abundances 
log fiNLTE = 5.28 + 0.11 dex and log s^^^ = 5.23+0.1 dex. The 
standard deviations cr are quoted as eiTors. 

It is important to note that there is no single set of measured 
oscillator strengths for all lines in our list. Data from five sources 
were used, and this undoubtedly reduces the precision of our re- 
sults, leading to a large standard deviation of the mean abun- 
dance. As a test, we calculated the mean log ej^^™ for each set 
of ^/-values (Table|2]i. The result is a large spread of abundances 
from 5.23 to 5.46 dex, which is due to strongly discriminate os- 
cillator strengths for the same lines measured by different au- 
thors. Indeed, for the majority of lines in our list the mean dif- 
ference log ^/(Oxford - Kurucz) = 0.096 dex, log ^/(Oxford - 
Becker) - 0.09 dex. There is, surprisingly, much better agree- 
ment between the absolute oscillator strengths of Becker et 
al. (119801) and Kurucz (1988 ): the mean difference (Oxford - 
Becker) - 0.047 dex. We have, nevertheless, serious doubts 
about the accuracy of all measurements: the errors, quoted by 
the authors (0.04 to 0. 1 dex), are much smaller than the actual 
differences between the ^/-values from these sources. 

Only 14 lines in our list (marked in the Table [1] with an as- 
terisk) have consistent measured oscillator strengths. When only 
those lines are used, we obtain logej^"^™ - 5.30 + 0.09 dex. 
When two obvious "outsiders" with log s =5.15 dex from mul- 
tiplet 16 are ignored, the abundance of manganese increases to 
5.36 + 0.06 dex. Despite its large standard deviation, this value 
based on only 12 lines is more reliable than the one calculated 
with all lines from our sample. 

As can be seen from Table [1] there is a large line-to-line 
spread of abundances. To investigate this we calculated mean 
abundances logs and standard deviations for each multiplet 
(Table [3]l. The abundances obtained with oscillator strengths 
from different sources are marked with an asterisk in the abun- 
dance column; these values are also weighted according to the 
errors in oscillator strengths. A plot of log e is shown in Fig. 
I^^. There is a large difference between log e for multiplets with 
E\o^f, = 2 . . . 3 eV and all other lines. It is interesting that a similar 
iiTegularity was noted by Blackwell et al. ( 1982 ) and Simmons & 
Blackwell (1 19821 1 from the analysis of Fe i lines in the solar spec- 
trum. Lines with excitation energies 2.18 ... 2.2 eV gave lower 
abundances than lines with smaller or larger isiow Inaccurate 
damping parameters (Unsold formula) or a vague dependence 
on the multiplicity of a term were suggested as the cause. 



^ All mean abundances in this paper are weighted according to the 
uncertainties in the g/-values for all lines; higher weights are assigned 
to lines with smaller errors 
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Fig. 4. Abundances of the 39 Mn i lines in Tabled (a): LTE (cir- 
cles) and NLTE (filled triangles) abundances. The mean LTE 
abundance is represented by dots, NLTE by a continuous line, 
(b): Mean abundances in multiplets. Filled circles refer to the 
standard NLTE calculations. Open triangles refer to the NLTE 
calculations with adjusted logC6 values (see text). The solid 
black horizontal line indicates the mean multiplet NLTE abun- 
dance, logfimean - 5.28, for all Hncs without logCe correction, 
the dotted horizontal line represents a corresponding value of 
loge,nean = 5.31 dcx after log Ce adjustment. Multiplet numbers 
are indicated, (c): NLTE abundances as a function of oscillator 
strengths. Filled circles and a continuous linear regression line 
denote Oxford oscillator strengths of Booth et al. ( |1984al l, open 
triangles and a dotted line refer to Kurucz's database. 
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Table 3. NLTE abundances log e of Mn i multiplets. Values ob- 
tained with log ^/-values from different sources are marked with 
an asterisk; A logs values are given relative to the weighted 
mean NLTE abundance log e^^^ - 5.28 dex. 



Mult. 




^low 


logCe 


l0g£ 


cr 


Alog£ 






[eV] 




[dex] 






1 


2 


0.0 


-31.4 


5.35 





0.07 


4 


5 


2.14 


-31.41 


5.42* 


0.04 


0.14 


5 


3 


2.15 


-31.0 


5.38 


0.06 


0.1 


16 


3 


2.3 


-30.8 


5.20 


0.09 


-0.08 


20 


1 


2.91 


-30.83 


5.35 




0.07 


21 


7 


2.9 


-30.74 


5.21* 


0.09 


-0.07 


22 


7 


2.92 


-30.67 


5.21 


0.10 


-0.07 


23 


1 


2.94 


-30.4 


5.58 




0.2 


27 


3 


3.06 


-30.64 


5.25* 


0.08 


-0.02 


28 


2 


3.07 


-30.7 


5.40* 


0.13 


0.12 


32 


2 


3.12 


-30.73 


5.31* 


0.11 


0.03 


49 


3 


4.42 


-31.6 


5.40 


0.02 


0.12 



The only common feature between our analysis and the iron 
results is the use of oscillator strengths measured by the Oxford 
group. These laboratory data cover 61% of all lines in our list; 
in particular, abundances for multiplets 16 and 22 are calculated 
solely with \og gf values from Booth et al. ( |1984al ). Standard 
deviations for these multiplet abundances are no larger than in 
other multiplets, although Fig. |4j; shows a weak correlation of 
abundances with the oscillator strengths, where the correlation 
coefficient is equal to -0.7. On this plot, abundances calculated 
using Booth et al.'s (I1984al ) data are compared with those calcu- 
lated with Kurucz's data. The linear regression curves are shown 
for both sets. However, the observed behaviour could be a ran- 
dom variation misinterpreted for a trend, or it may be the conse- 
quence of incorrectly chosen damping parameters. 

In order to test the influence of the Van der Waals damp- 
ing we carried out spectrum synthesis calculations adjusting Ce 
for each multiplet. The correction of this parameter was found 
from the requirement of obtaining equal abundances for weak 
and strong lines of a common multiplet; A log C(, for all inves- 
tigated multiplets averaged to -0.7 dex relative to the Anstee & 
O'Mara data, which is close to values predicted by the Unsold 
formula. However, only the lines of intermediate multiplet num- 
bers demonstrated a weak sensitivity to the Van der Waals damp- 
ing (Fig.ih): Aloge(16)= -1-0.14 dex, Aloge(21)= -h0.03 dex, 
Aloge(22)= -1-0.04 dex, Aloge(27)= -1-0.03 dex. This is insuf- 
ficient to explain the discrepancy of 0.2 dex between the abun- 
dances from lines of these multiplets and the others. Besides, it 
is not always possible to discern the influence of Ce and log e 
on the line profiles. This is demonstrated in Fig. |5] where three 
combinations of both parameters, i.e., smaller abundance and 
larger damping, or larger abundance and smaller damping, lead 
to equally good fits. The difference between the profiles is visi- 
ble only in the blend on the blue line wing. In our final estimate 
of the solar Mn abundance we decided to keep the Anstee & 
O'Mara damping parameters, even though they require a signif- 
icantly smaller abundance for particular lines. 

The discrepancy between the abundances (both NLTE and 
LTE) from different multiplets is hard to explain. We do not be- 
lieve that errors in damping could be as large as this, neither 
with confidence can we attribute this to errors in the oscilla- 
tor strengths. It is interesting that a similar result was obtained 
for Mm lines by the other authors. Lines that give a strong un- 
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Fig. 5. NLTE line profiles of the line 4823 A. Different combina- 
tions of damping parameters and resulting abundance deviations 
from log SMn.G = 5.47 dex are shown. 



1.0 - 




Fig. 6. LTE line profiles of the Mn ii line at 3488 A. The sti'ong 
line wings centered near 3489.05 A are most likely due to an 
autoionization line, represented as an Fei transition at 3.5 eV 
with an extremely large radiative damping constant. 



derabundance in our analysis also have smaller log e relative to 
the average value in the work of Becker et al. (1980). Booth et 
al. ( 1984b) also distinguished two groups of lines based on the 
lower excitation potential: the eV lines, which show consider- 
ably higher abundance, and those between 2 - 3 eV with lower 
abundances, but are sensitive to damping enhancement. To re- 
solve this inconsistency we need new quality measurements of 
oscillator strengths and, perhaps, a depth-dependent microturbu- 
lence or a 3D hydrodynamical model. 

Unfortunately, we can not solve the main problem: there is a 
strong inconsistency between the solar photospheric abundance 
of Mn and a meteoritic value of 5.50 ± 0.03 dex (Lodders 2003I I. 
Quite on the contrary to the expectations, every subsequent anal- 
ysis with more precise atomic data and model assumptions gives 
progressively decreasing log SMn,©- Generally, high values that 
are consistent with the meteoritic abundance (e.g., Greenlee & 
Whaling |1979t Becker et al. 119801 ) result from erroneous os- 
cillator strengths, neglect of HFS, an LTE approach (our work 
is t\\t first NLTE analysis), and/or unaccounted blends. For in- 
stance, Greenlee & Whaling (119791 ) used only six lines with 
Wa < 14 mA, assuming that there are no unrecognized blends 
in the solar spectrum and neglecting HFS. Yet, we have found 
that four of these six lines contain blends, we therefore ignored 
them for our analysis. The most recent and commonly used 
LTE abundance of Mn (5.39 + 0.03 dex) is that of Booth et al. 
C 1984b) . It is also rather low, hence the authors suggested that 
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NLTE effects could be the cause of methodical errors. Indeed, 
NLTE effects of Mn i in the Sun are significant, but they only in- 
crease the existing discrepancy. Even if only the weak lines with 
Wa < 50 mA are used, which are relatively insensitive to the 
keystones of the abundance analysis (damping, model structure, 
HFS, microturbulence), we obtain a weighted NLTE abundance 
of 5.37 ± 0.06 dex. The incorrect abundance of the reference 
element Si may be the cause. Some authors compare the abun- 
dance of manganese with that of iron. The recent NLTE esti- 
mate is logSFe,0 = 7.48 . . .7.51 dex (Gehren et al. 1200 IbK If 
we adopt this value, a ratio Fe/Mn = 148 is derived; it is sig- 
nificantly larger than Fe/Mn = 89 in CI chondrites (Cameron 
|1973t . Perhaps, as suggested by Booth et al. ( 1984b ), the chem- 
ical composition of CI meteorites does not fully represent the 
composition of the proto-Sun. 

To investigate the ionization equilibrium of manganese, 
we have tried to determine the abundance from Mnn lines. 
Martinson et al. ( 119771 ) claim that Mn ii lines can be used to ob- 
tain the same photospheric solar abundance as found from Mn i 
lines. From equivalent width measurements of 9 lines of ionized 
manganese these authors have calculated logfiMn.o = 5.4 + 0.2 
dex. We have checked the solar spectrum of Mnii and found 
that no lines are suitable for any of the abundance determination 
methods. All lines strong enough to be discerned in the solar 
spectrum are located in the near-UV, and are heavily blended, 
which complicates the continuum placement. We have, never- 
theless, synthesized these lines under the assumption of LTE: 
all lines of interest have saturated cores, i.e. the variation of 
the abundance within + 0.2 dex affects only the strength of the 
wings. Thus, due to the presence of strong blends in the line 
wings, no conclusion can be made about the true abundance. For 
the line at A 3488 A the effect of HFS was investigated (see Fig. 
|6l). Data for the levels involved in this transition were taken from 
Holt et al. (1 19991 ). Hyperfine structure broadens the profile and 
leads to a better fit of the inner wings provided the microturbu- 
lence is reduced to - 0.8 km s"'. Its effect is very significant: 
a line profile calculated with no HFS and log s^^^ - 5.44 dex is 

equal to a profile with 2 HFS components and \ogE^^^ - 5.29 
dex. Hence, abundances for the Mnii lines calculated without 
accounting for HFS can easily be overestimated by 0.2 dex. Our 
results bring into perspective the conclusions of Martinson et al. 
( 119771 1 : the neglect of HFS in their calculations is, most likely, 
responsible for an overestimated abundance of manganese. At 
the same time, the comparison of abundances derived from two 
ionization stages of manganese turns out to be an unreliable tech- 
nique. It should be performed with caution. 
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Fig. 7. Synthetic NLTE profiles for the Mm lines, (a): Variation 
of /15420 A with photoionization cross-sections multiplied by 
different scaling factors 5p. (b): Response of /16016 A to vari- 
ous scaling factors for Hi collisions, (c): The Mm line 8703 A 
displaying the dependence on the HFS constant A{ J). The line 
suffers from an unknown blend. 



4.3. Uncertain parameters 

4.3.1 . Photoionization and collisions with H i atoms 

We performed NLTE calculations with various scaling factors 
for photoionization and collision cross-sections. The effect on 
the departure coefficients is described in detail in Sect. 3.2. Here 
we show how changing these parameters affects the line profiles 
and the resulting abundances. 

The behaviour of profiles for different scaling factors Sp to 
hydrogenic photoionization cross-sections is shown in Fig. |7^. 
The difference in the profiles is only marginal: maximum abun- 
dance corrections amount to -0.04 dex with Sp increasing from 
to 1000. This effect was confirmed with twelve lines. We con- 
clude that the photoionization cross-sections are not crucial in 
abundance calculations as long as very large Sp enhancements 



(by a few orders of magnitude) to hydrogenic cross-sections are 
avoided. 

A similar test was performed for collisions with H i atoms, 
the parameter usually considered the most uncertain in NLTE 
calculations. Scaling factors Su - 0.05, 0.5, 1,5 and 10 were ap- 
plied to b-b and b-f cross-sections computed from Drawin's for- 
mula. The NLTE profiles of the line at 6016 A are presented for 
five different assumptions in Fig.|7]5- As expected, with increas- 
ing rates of collisional excitation and ionization by H i atoms, 
the NLTE profile approaches LTE. However, the actual coiTec- 
tion to the abundance is rather small: for ten lines investigated it 
does not exceed 0.06 dex for 5h = 5. Neglect of hydrogen colli- 
sions (not shown here) requires a reduction of the abundance by 
0.01 dex. According to the results of Belyaev & Barklem (20031, 
Drawin's formula (S-a - 1) strongly overestimates the rate of 
transitions due to collisions with H i. Therefore, this scaling fac- 
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Table 4. Average depths of Mn i NLTE line formation for a line 
with continuum according to Achmad et al. ( |1991| l. 
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tor is already too large, and we conclude that the magnitude of 
uncertainty concerned with the choice of Sh is not much larger 
than 0.05 dex. 



4.3.2. Hyperfine splitting data 

The sensitivity of lines to uncertainties in the hyperfine split- 
ting data was investigated (Fig. [T};). It is seen that differences 
in the 8703 A line profiles computed for A{J) = 0.015 cm"', 
A(J) = 0.0158 cm ', and A(7) = 0.0175 cm ' are rather small. 
These differences correspond to an abundance correction of 0.02 
dex. A notable change in the profile is found only when A{J) is 
modified by more than 4 mK (ImK = 0.001 cm"'). However, 
the precision of A{J) measurements that we use is very high; the 
average error is 0.5 mK. We have tested other lines, including 
those with saturated cores, varying the A{J) value for the cor- 
responding levels by ~ 2 mK. In this case, the relative shift of 
the components leads to A log s ~ 0.01 - 0.02 dex. We conclude 
that due to the much stronger influence of the other parameters 
(log gf, Cg), uncertainties in HFS data do not affect the resulting 
abundances. However, correcting the A(J) value may in a few 
cases help to obtain a better fit to the observed spectral lines. 

The influence of HFS on line formation depths is interesting. 
The eff'ect of including HFS on the position and width of forma- 
tion layer for the weak resonance line 5394 A was investigated 
by Vitas ( |2005l l. The main result was that the line was formed 
higher in the photosphere when HFS was not taken into account. 
We have tested this observation on a sample of lines with differ- 
ent Wa and gf values; their depths of formation computed with 
and without HFS are given in Table |4] These are average depths 
of NLTE line formation calculated with the contribution func- 
tion to the total emergent radiation according to Achmad et al. 
( 1199 1). Both for strong and weak lines the depth of core for- 
mation decreased when hyperfine splitting was neglected in line 
synthesis. The corresponding shift in local temperatures Ar is 
given in the last column. We note that the formation of weaker 
lines takes place at significantly higher Tg. For instance, the res- 
onance line at 5394.67 A with 6 HFS components is formed 
in the layer with = 5340 K, rather than in the layer with 
Te - 4920 K (depth of formation for the same line with a single 
component). This result shows that weak lines are to a greater 
degree affected by inclusion of HFS in line synthesis, compared 
to stronger lines. 



5. Conclusions 

The NLTE formation of Mm lines in the solar atmosphere is 
rather complex. Our results indicate that manganese belongs to 
a class of mixed domination atoms, where several channels com- 
pete in establishing non-equilibrium level populations. However, 
in our model neither of the processes can be regarded as a dom- 
inating type: in particular, the absence of photoionizations does 
not restore an LTE distribution of populations. We suggest that 
the photoionization inefficiency is the simple result of too low 
cross-section approximations. Whether with or without quantum 
defect, the hydrogenic approximations deliver values roughly 
factors of 100 to 1000 lower than seen in comparable atoms for 
which better calculations are available (Fe, Si). 

The depopulation of the Mm ground state and low excita- 
tion levels in the absence of photoionization can be understood 
as follows: radiative transitions with strong optical pumping ef- 
fects cause an overpopulation of the high-excitation levels with 
respect to the low-excitation levels. Due to the strong coUisional 
coupling of the high-excitation levels with the Mn ii ground 
state, which is in LTE, the departure coefficients of these lev- 
els will be close to unity; consequently the departure coefficients 
of the low-excitation levels will drop significantly below unity. 
In the absence of the high-excitation levels. Mm behaves as a 
photoionization-type ion. Disregarding inelastic collisions with 
electrons results in a weaker coupling of the levels to each other 
and a stronger deviation of line source functions from By{Ts). 

All levels that are involved in transitions of interest are un- 
derpopulated relative to LTE. Hence the NLTE assumption gen- 
erally leads to an increase of the Mn i depth of line formation 
that forces a weakening of the lines. However, for every partic- 
ular line the ratio of the upper and lower level population, de- 
termined by an interplay of the above-defined NLTE processes, 
modulates the behaviour of the wings and cores, and hence the 
value and sign of NLTE corrections. For the weak and interme- 
diate strength lines (Wa < 80 mA), NLTE and LTE assumptions 
give a similarly good fit to the observations, provided a certain 
abundance correction is performed. The NLTE corrections for 
these lines are positive and amount to 0. 1 dex. Formation of the 
stronger lines (Wa > 80 mA) is characterized by an amplified 
absorption in the core and a decreased absorption in the wings, 
compared to LTE. NLTE corrections for strong lines are scat- 
tered around zero or negative. 

The LTE analysis of Mm has two pitfalls: the strong lines 
can not be reproduced due to different profile shapes, and abun- 
dances of the weak lines under LTE are overestimated. For this 
reason we do not recommend the commonly used solar LTE 
abundance logeMn.o = 5.39 + 0.03 dex, obtained by Booth et 
al. ( 1984b). Instead, an NLTE abundance of 5.36 + 0.1 dex, de- 
termined from 12 lines of Wa = 12-140 mA with relatively 
reliable oscillator strengths in the current analysis, is suggested. 
We consider abundances from other lines in our list as unreliable 
due to (1) high sensitivity of particular lines to Ce (for multiplet 
16 Alogs = H-0.14 dex for AlogCg = -0.9) and (2) incon- 
sistent oscillator strengths between different sources. This result 
increases the discordance with the abundance of Mn in CI mete- 
orites of 5.50 + 0.03 dex (Lodders |2003] l. 

The atomic parameters influence neither the line profiles nor 
the derived abundances significantly. The assumption of hydro- 
genic photoionization cross-sections leads to an overestimation 
of the abundance by 0.04 dex in comparison to the case when the 
enhancement factor of 1000 is applied. Despite the importance 
of the HFS, no significant errors associated with the choice of 
HFS data are expected in the analysis. The 0.5 mK accuracy of 
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measurements produces an uncertainty of 0.01 dex in the abun- 
dance. Anstee & O'Mara's damping constants for stronger lines 
in some multiplets lead to abundances that are smaller than those 
obtained from the weaker lines. Adjustment of C(, values inside 
each multiplet increases the average abundance by 0.03 dex. 
However, the discrepancy of 0.2 dex between the abundances 
from lines of multiplets 16, 21 and 22 and all other multiplets 
still remains. We can not with confidence attribute this to errors 
in damping or in the oscillator strengths. Two findings provide 
evidence for the latter: (1) correlation of abundances with the 
Oxford oscillator strengths and (2) similar regularity found by 
Simmons & Blackwell et al. ( fT9g2t and Blackwell et al. (1982*) 
in the analysis of solar Fei lines with Oxford log gf values. 
Collisions with neutral hydrogen atoms are able to thermalize 
the levels only if an enhancement factor > 10, with respect 
to Drawin's cross-sections, is applied. For the smaller enhance- 
ment factors, although the differences between NLTE and LTE 
profiles decrease, the abundance corrections do not exceed 0.05 
dex. Total neglect of collisions with hydrogen requires a reduc- 
tion of the abundance by 0.01 dex. 

The main source of errors in this work is the uncertainties 
of measured oscillator strengths, which are inaccurate within a 
range of 0.04 to 0.1 dex. However, only the product of abun- 
dance and oscillator strengths log gfs will be used for studies of 
manganese in metal-poor stars, where deviations from LTE are 
expected to be large. 
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Fig. 8. Grotrian diagram of the Mn i model atom. Solid lines represent allowed and forbidden transitions included in the model 
atom. 
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Table 5. HFS constants A and B (in units of 10 cm ') representing magnetic dipole and electric quadrupole interactions for Mm 
levels. Level energies E are given in eV. 
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